function DB_askari = generate_DB_askari()
    % Correlations for pure plasmas
    %
    % Coefficients to calculate frozen and equilibrium specific heat and
    % constant pressure, specific enthalpy, specific entropy and mean molar
    % mass according to [1]
    %
    % [1] Askari, O., 2017, “Thermodynamic Properties of Pure and Mixed
    %     Thermal Plasmas Over a Wide Range of Temperature and Pressure,”
    %     ASME J. Energy Resour. Technol., 140(3), p. 032202.
    %
    % Species included: H2, He, C, N2, O2, Ne, Ar
    %
    % Returns:
    %     DB_askari (struct): Database

    DB_askari.H2 = generate_H2_plasma();
end

% SUB-PASS FUNCTIONS
function species = generate_H2_plasma()
    % Correlation for H2 pure plasma

    % I = [I1, I2, I3]
    species.I = [2, 2, 2];
    % m = [m]
    species.m = 5;
    % coeffi = [zeta0_i, zeta1_i, zeta2_i, zeta3_i, zeta4_i, zeta5_i]
    species.a_f = [1.884e+00, 2.549e-02, 9.699e-04, -7.532e-04, -1.115e-04, -4.259e-06; ...
                3.008e+00, -1.453e-03, 7.572e-05, 4.345e-05, 5.105e-06, 1.854e-07];
    species.b_f = [7.802e+00, 1.931e-02, -1.205e-03, -2.162e-04, -1.420e-05, -2.728e-07; ...
                    1.509e-06, 3.586e-05, 1.011e-04, -1.769e-03, 4.694e-02, 9.573e+00];
    species.c_f = [3.730e-06, 9.551e-05, 7.664e-04, 3.146e-03, -6.577e-02, 1.088e+00; ...
                    -1.554e-05, -4.362e-04, -3.287e-03, 2.152e-03, 4.540e-02, 2.541e+00];
    species.a_r = [-9.914e-07, -2.866e-05, -3.110e-04, -4.130e-03, -1.313e-01, 5.075e+00; ...
                    6.298e-07, 1.186e-05, -8.989e-05, -4.773e-03, -1.365e-01, 5.491e+00];
    species.b_r = [-3.446e-08, 5.403e-07, 6.375e-05, 2.005e-03, 6.572e-02, 8.238e+00; ...
                    2.037e-07, 6.145e-06, 8.638e-05, 1.970e-03, 7.367e-02, 9.617e+00];
    species.c_r = [1.029e-06, 2.974e-05, 3.031e-04, 2.717e-03, 6.785e-02, -1.535e+00; ...
                    3.729e-07, 1.301e-05, 1.899e-04, 2.297e-03, 6.470e-02, -1.420e+00];
    species.a_m = [-5.231e-08, -1.738e-06, -2.123e-05, -1.294e-04, -4.359e-04, 1.042e-02; ...
                    3.892e-08, 1.223e-06, 1.536e-05, 1.216e-04, 5.973e-04, -6.903e-01];
    species.b_m = [1.450e-07, 5.093e-06, 9.390e-05, 2.003e-03, 6.444e-02, 8.191e+00; ...
                    -8.348e-08, -8.065e-07, 5.746e-05, 2.254e-03, 7.428e-02, 9.565e+00];
    species.c_m = [1.143e-07, 1.780e-06, -4.486e-05, -1.977e-03, -6.380e-02, 2.355e+00; ...
                    2.053e-07, 2.706e-06, -8.753e-05, -2.744e-03, -6.868e-02, 2.229e+00];
    species.lambda = [1.313e-05, 3.261e-04, 1.949e-03, -4.339e-03, -6.728e-02, 1.456e+01; ...
                    1.015e+00, 2.672e+01, 1.796e+02, -1.131e+02, 3.337e+03, 9.082e+05; ...
                    1.336e-04, 3.692e-03, 2.945e-02, 1.008e-02, -1.151e+01, -1.808e+02; ...
                    0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 2.016e+00];
end

% species.a1_f = [1.884e+00, 2.549e-02, 9.699e-04, -7.532e-04, -1.115e-04, -4.259e-06];
% species.a2_f = [3.008e+00, -1.453e-03, 7.572e-05, 4.345e-05, 5.105e-06, 1.854e-07];
% species.b1_f = [7.802e+00, 1.931e-02, -1.205e-03, -2.162e-04, -1.420e-05, -2.728e-07];
% species.b2_f = [1.509e-06, 3.586e-05, 1.011e-04, -1.769e-03, 4.694e-02, 9.573e+00];
% species.c1_f = [3.730e-06, 9.551e-05, 7.664e-04, 3.146e-03, -6.577e-02, 1.088e+00];
% species.c2_f = [-1.554e-05, -4.362e-04, -3.287e-03, 2.152e-03, 4.540e-02, 2.541e+00];
% species.a1_r = [-9.914e-07, -2.866e-05, -3.110e-04, -4.130e-03, -1.313e-01, 5.075e+00];
% species.a2_r = [6.298e-07, 1.186e-05, -8.989e-05, -4.773e-03, -1.365e-01, 5.491e+00];
% species.b1_r = [-3.446e-08, 5.403e-07, 6.375e-05, 2.005e-03, 6.572e-02, 8.238e+00];
% species.b2_r = [2.037e-07, 6.145e-06, 8.638e-05, 1.970e-03, 7.367e-02, 9.617e+00];
% species.c1_r = [1.029e-06, 2.974e-05, 3.031e-04, 2.717e-03, 6.785e-02, -1.535e+00];
% species.c2_r = [3.729e-07, 1.301e-05, 1.899e-04, 2.297e-03, 6.470e-02, -1.420e+00];
% species.a1_m = [-5.231e-08, -1.738e-06, -2.123e-05, -1.294e-04, -4.359e-04, 1.042e-02];
% species.a2_m = [3.892e-08, 1.223e-06, 1.536e-05, 1.216e-04, 5.973e-04, -6.903e-01];
% species.b1_m = [1.450e-07, 5.093e-06, 9.390e-05, 2.003e-03, 6.444e-02, 8.191e+00];
% species.b2_m = [-8.348e-08, -8.065e-07, 5.746e-05, 2.254e-03, 7.428e-02, 9.565e+00];
% species.c1_m = [1.143e-07, 1.780e-06, -4.486e-05, -1.977e-03, -6.380e-02, 2.355e+00];
% species.c2_m = [2.053e-07, 2.706e-06, -8.753e-05, -2.744e-03, -6.868e-02, 2.229e+00];
% species.lambda1 = [1.313e-05, 3.261e-04, 1.949e-03, -4.339e-03, -6.728e-02, 1.456e+01];
% species.lambda2 = [1.015e+00, 2.672e+01, 1.796e+02, -1.131e+02, 3.337e+03, 9.082e+05];
% species.lambda3 = [1.336e-04, 3.692e-03, 2.945e-02, 1.008e-02, -1.151e+01, -1.808e+02];
% species.lambda4 = [0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 0.000e+00, 2.016e+00];
